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ABSTRACT 

Through  the  use  of  a  predictive  model  and  experimentation  in 
a  calorimeter,  the  temperatures  and  humidities  which  represent 
equal  comfort  states  to  the  human  for  minimal  ventilation  rates  in 
a  neutral  thermal  radiation  environment  are  established.  This  is 
accomplished  by  establishing  from  experimental  data  a  relationship 
between  the  subject's  skin  temperature  and  the  per  cent  01  sweat 
glands  that  are  activated.  Points  on  the  resulting  curve  are 
defined  as  values  on  the  Physiological  Thermal  Index.  They  are 
used  in  conjunction  with  experimentally  obtained  heat  and  mass 
transfer  coefficients  to  place  lines  on  the  psychrometric  chart 
which  represent  constant  values  of  this  index.  The  effects  of 
variations  in  air  flow  rate,  and  metabolic  output,  are  presented. 
The  predicted  results  are  also  compared  with  experimental  evidence 
from  the  calorimeter  and  with  the  effective  temperature  values 
from  the  literature.  The  evidence  to  date  indicates  that  the 
predictive  model  is  sound  and  that  the  Physiological  Thermal  Index 
affords  a  new  and  clearer  way  to  view  the  interaction  between  man 
and  his  thermal  environment. 

A  predictive  model  for  determining  human  thermal  response 
to  transient  radiant  flux  of  the  type  associated  with  building 
fires  in  spaces  adjacent  to  survival  shelters  is  also  present. 

It  is  built  upon  the  concepts  previously  discussed  and  includes 
the  thermal  capacitance  of  the  test  subject.  Depending  upon  the 
onset  rate  of  flux  generation,  it  is  shown  that  the  skin  tem¬ 
perature  can  be  driven  above  body  core  temperature.  A  'plan  for 
experimentally  verifying  the  results  from  the  predictive  model 
using  the  human  calorimeter  is  discussed. 
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I.  INTRODUCTION 

Previous  Work  Summary 

This  investigation  was  undertaken  as  a  continuation  of  the 
Monoman  Calorimeter  Project  Cl]*,  which  is  a  study  of  the  human 
thermal  reaction  to  simulated  crowded  conditions  and  low  ventila¬ 
tion  rates.  For  the  past  five  years  tests  have  been  conducted  on 
a  sample  population  of  adult  males ,  females  and  children  under 
high  dry  bulb  temperatures,  high  humidities  and  very  low  air  flow 
rates.  Most  of  the  data  generated  was  along  82,  85,  and  90° F  effec¬ 
tive  temperature  lines.  The  previously  reported  work  primarily 
examined  metabolic  rates  and  the  process  lines  generated  by  test 
subjects  in  this  severe  thermal  environment.  Based  on  the  experi¬ 
mental  evidence  it  was  shown  that  a  standard  population  sample 
under  sedentary  conditions  generates  an  average  metabolic  output 
of  275  Btu/hr.  This  indicated  that  the  400  Btu/hr  commonly  jsed 
in  planning  has  a  considerable  safety  margin  unless  the  shelterees 
are  all  young  men  who  indeed  do  have  a  typical  metabolic  rate  of 
400  Btu/hr.  The  concepts  of  non-isostate  and  isostate  environ¬ 
ments  were  established  as  the  extreme  process  limits  between  which 
an  actual  shelter  process  must  lie. 

During  this  earlier  work  part  of  the  investigation  was  to 
evaluate  the  effective  temperature  lines,  developed  by  F.  C. 

Houghton  and  C.  P.  Yaglou  [2],  from  subjective  response  of  the 
test  subjects.  The  question  being  whether  the  effective  tem¬ 


perature  concept  would  truly  represent  equal  comfort  states  for 
minimal  air  circulation.  Many  laboratories  and  field  workers  [3] 


The  numbers  m  brackets  designate  references  listed  in 
the  bibliography. 
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have  found  that  this  index  overestimates  the  influence  of  humid¬ 
ity  on  sensations  of  warmth  and  comfort  at  ordinary  temperatures, 
and  underestimates  the  effect  in  very  high  temperatures.  The 
effectKe  temperature  lines  indicate  a  person’s  feeling  of 
warmth  immediately  after  entering  a  conditioned  space.  Over¬ 
estimation  of  "he  humidity  effect  in  the  absence  of  sweat  gland 
activity  (referred  to  here  as  insensible  evaporation)  is  explain¬ 
ed  by  adsorpt Ion  and  desorption  phenomena  and  by  failure  to  take 
adaptation  into  account  [4].  Several  studies  have  been  done  to 
improve  the  ET  lines.  Hardy  and  DeBois  [5],  and  Gagge,  et  al.  , 
[6]  have  found  a  close  relationship  between  mean  skin  temperature 
and  comfort  in  air  conditions  that  lie  below  the  zone  of  active 
evaporative  cooling.  C.  P.  Yaglou  has  found  mean  skin  tempera¬ 
ture  is  a  fairly  sensitive  index  of  warmth  in  temperatures  under 
82°F,  or  when  persons  are  not  actively  sweating.  He  suggested 
that  the  mean  skin  temperature  be  taken  as  an  index  of  comfort 
when  there  is  no  active  evaporative  cooling.  Tn  the  zone  of 
active  evaporation  cooling  the  skin  temperature  rises  very  slowly 
w.iile  the  evaporative  heat  loss  rises  very  sharply  [73. 

Many  experiments  have  been  made  to  study  the  relation 
between  total  evaporative  heat  loss  and  operative  temperature 
when  radiant,  convective  and  evaporative  heat  exchanges  are 
present.  C.  E.  Winslow,  L.  P.  Herrington  and  A.  P.  Gagge 
studied  the  relation  between  total  evaporative  heat  loss  and 
mean  skin  temperature  [7]. 

Present  Program  Objectives 

(a)  To  establish  a  physical  basis  for  the  prediction  of  lines 
of  equal  comfort  for  semi-nude  persons  in  equilibrium  with 
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environmental  conditions  of  very  low  air  flow  rates  and  with¬ 
in  the  range  wherein  the  human  body  can  maintain  a  thermal 
regulatory  control. 

(b)  To  use  this  as  a  basis  for  predicting  lines  of  equal  comfort 
on  the  psychrometric  chart. 

(c)  To  compare  test  subject  response  with  the  predictive  model  by 
use  of  the  calorimeter. 

(d)  To  extend  the  predictive  model  to  include  transient  behavior 
associated  with  a  thermal  radiation  pulse. 

Physical  Basis  for  the^Predictive  Model 

The  physical  basis  for  the  predictive  model  lies  in  three 
relationships: 

(a)  The  relationship  between  the  skin  temperature  and  percentage 
of  the  sweat  glands  that  are  actively  sweating. 

(b)  Energy  balance  equations  between  the  human  body  and  the 
environment. 

(c)  Heat  and  mass  transport  rate  equations  between  the  person  and 
his  surroundings. 

Experimental  Apparatus 

The  apparatus  used  to  gather  the  experimental  data  is  the 
monoman  calorimeter,  shown  in  Fig.  1,  which  was  designed  original¬ 
ly  to  simulate  the  environment  of  an  overcrowded  survival  shelter 
under  the  condition  of  very  low  air  flow  rates.  The  main  compo¬ 
nents  of  the  calorimeter  are  the  adiabatic  test  cell  in  which  the 
test  subject  is  encapsulated  and  the  ventilation  system.  The  test 
cell  is  a  fabricated  bag  made  of  polyethylene  sheet  and  plastic 
tape.  The  bag  is  lined  on  the  inside  with  aluminized  mylar  while 


MONOMAN  CALORIMETER 


FIGURE  i 


it  is  covered  on  the  outside  with  several  sheets  of  aluminized 
mylar  separated  by  a  sheet  of  air  cap.  The  net  radiative  heat 
exchange  between  the  calorimeter  and  the  test  subject  is  negli 
ble  because  of  the  reflective  walls  of  the  calorimeter  which 
simulates  being  surrounded  by  others.  The  bag  has  excellent 
thermal  response  characteristics  and  high  thermal  resistance. 
Adiabatic  conditions  are  maintained  by  automatically  adjusting 
the  intensity  of  five  rings  of  lamps  that  radiate  heat  to  the 
calorimeter  wall.  Air  supply  for  the  calorimeter  with  desired 
quantity,  temperature  and  humidity  is  from  an  air-conditioning 
apparatus.  For  more  details  see  reference  Cl]. 

Instrumentation  and  Test  Procedure 


Measured  quantities  of  primary  interest  are  wet  and  dry  bulb 
temperatures  at  inlet  and  exhaust,  air  flow  rate,  and  temperature 
gradients  across  the  walls  and  base  to  insure  adiabatic  condition. 
Thermocouples  are  used  to  measure  wet  and  dry  bulb  temperatures 
while  a  positive  displacement  flow  meter  and  stop  watch  are  used 
to  measure  the  air  flow  rate. 

Physiological  measurements  of  primary  interest  are  pulse, 
oral  temperature,  and  skin  temperature.  The  skin  temperature  is 
taken  as  the  weighted  average  of  readings  of  nine  thermistors 
located  at  different  points  on  the  body  surface.  As  the  skin 
temperature  harness  was  developed  as  part  of  the  c -‘'rent  project, 
it  is  discussed  in  detail  in  Appendix  1. 

After  the  system  is  brought  to  the  desired  condition,  the 
person  to  be  tested  is  secured  in  the  calorimeter.  The  equilib¬ 
rium  condition  of  the  person  and  the  calorimeter  is  reached  by 
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automatically  adjusting  the  heating  lamps  until  the  calorimeter 


walls  are  adiabatic  while  holding  air  inlet  temperature,  humidity 


and  flow  rate  constant.  Data  is  taken  after  observing  constant 


conditions  of  the  calorimeter  and  the  subject. 


An  analytical  attempt  to  correlate  the  convective  coeffi¬ 


cients  with  low  air  velocities  and  with  the  free  convection  para¬ 


meters  [8]  failed  because  the  random  motions  of  the  test  subject 


dominate  over  the  fluid  motion  produced  by  the  free  convection 


mechanisms.  It  is  for  this  reason  that  the  experimentally  deter¬ 


mined  convective  heat  and  mass  transfer  coefficients  from  many 


calorimeter  experiments  are  averaged  and  called  natural  motion 


limit  convective  heat  and  mass  transfer  coefficients. 


Report  Presentation  Format 


Section  II  will  deal  with  the  energy  balance  and  transport 


rate  equations.  From  these  equations  and  the  experimental  data 


generated  from  test  subjects,  average  values  of  convective  heat 


and  mass  transfer  coefficients  are  determined.  In  Section  III, 


.  ional  experimentally  determined  physiological  data  will  be 


examined  to  establish  the  criteria  of  human  comfort  which  form 


the  basis  for  prediction  of  equal  comfort  states  under  different 


environmental  conditions.  Section  IV  will  show  how  to  generate 


equal  comfort  lines  on  the  psychrometric  chart  for  various  condi¬ 


tions  of  air  flow  rate,  metabolic  rate,  temperature,  and  humidity 
based  on  constant  skin  temperature  and  constant  percentage  of  the 


body's  sweat  glands  that  are  active.  This  section  will  also 


examine  experimental  evidence  obtained  from  the  calorimeter  and 


other  experiments  in  the  light  of  these  predictions.  Section  V 
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will  show  the  extension  of  the  predictive  technique  to  include  the 
transient  thermal  behavior  of  a  shelteree  when  exposed  to  a  time 
dependent  radiation  flux.  The  energy  Storage  effect  of  the 
subject's  body  and  peripheral  tissues  and  the  thermal  resistance 
to  heat  flow  within  the  body  are  taken  into  account.  A  summary 
and  conclusions  relative  to  the  work  to  date  and  recommendations 
for  follow-on  work  constitute  Section  VI. 
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II.  ENERGY  BALANCE  AND  TRANSPORT  RATE  EQUATIONS 

1  .  ;  •  1 

Heat  and  Mass  Transfer  Rate,  Equations  .  1 

The  heat  balance  equation  describing  the  heat  exchange 
between  the  human  body  and  its'  surrounding?  is  'rased  on  the  law 

:  i  •  ■ 

of  conservation  of  energy.  This  can  be  stated  as:  the  net  rate 

i  i 

of  energy  flow  across  the  skin  plus  the  external  mechanical  work 

‘  I  X  \  l 

delivered  by  the  human  body  equals  the  metabolic  rate  minus  the 

•  i  * 

increase  in  heat  content  of  the  body.  This  basic  relation 

I 

applies  regardless  of  the  source  of  heat  or  the:  nature  of  its  . 


ii 
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transfer  [9].  This  statement  can, be  mathematically  expressed  as: 


MR-S  =  E  +  R+  C  +  W 


(1) 


E  = 
R  = 
C  = 
W  = 


where  MR  =  (Btu/hr)  the  net  rate  of  metabolic  heat  product ' oh 
S  =  (Btu/hr)  the  energy  storage  rate  of  the  body 

! 

(Btu/hr)  rate  of  evaporative  energy  loss . 

(Btu/hr)  rate  of  thermal  energy  loss  by  radiation 

1  • 

(Btu/hr)  rate  of  thermal  energy  loss  by  convection 

i 

(Btu/hr)  mechanical  work  accomplished  ! 

The  radiative  heat  output  is  negligible  in  the  monoman 

calorimeter  and  will  be  omitted  in  the  initial  stages  of  ,this 

«  : 

development.  The  energy  storage  rate  of  the  body1 is  taken  as 

j 

zero  initially  for  equilibrium  conditions  are  assumed  at  the  out-' 
set.  Therefore,-  the  energy  balance  equation  for  sedentary  sub¬ 
jects  (not  delivering  external  work)  can  be  reduced  to: 


*  1 


-23 

.^3 

1 


.a* 
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MR  =  C  +  E 


(2) 
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The  convective  heat  transport  rate  equation  can  be  expressed 


1  f 

V  1 


(3) 


i  ^ 


.e 
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where 

hc  =  (Btu/hr  ft^°F)  is  the  body  area  averaged  coefficient 
of  heat  transfer. 

A  =  (ft  )  is  the  human  body  surface  area  calculated  from 

DuBois  empirical  equation  [10]. 

T  =  (op)  -the  body  area  averaged  skin  temperature. 

^  S 

(See  Appendix  1  for  experimental  method  of  determining 
this  value.) 

T  =  (°F)  is  the  logarithmic  mean  average  of  air  inlet  and 

cl 

exhaust  temperatures  and  is  determined  by  the  follow¬ 
ing  equation  (see  Appendix  2): 


where 

Tq  =  Dry  bulb  temperature  at  exit  state  of  the  calorimeter. 

=  Dry  bulb  temperature  at  inlet  state. 

The  logarithmic  mean  average  is  used  in  calculating  T_ 
because  there  are  large  air  temperature  differences  between  inlet 
and  exhaust  states  when  test  subjects  are  exposed  to  inlet  dry 
bulb  temperatures  considerably  above  or  below  their  skin  tem¬ 
peratures  and  with  very  low  air  flow  rates. 

The  evaporative  energy  transfer  rate  equation  is 

E  =  hm  hfg  A  x  <ws  -  V  (5) 

o 

hm  =  (lb  water  vapor/hr  ft  (lb  water  vapor/lb  dry  air)) 

average  convective  mass  transfer  coefficient  based  on 
a  wet  body  area  of  (xA). 


where 


»  msx^StS^1. 


ft  I  i- 

it  i 

ii 

S  I  T 

if  I 

i  I 

1 1  t 

i  i 

ij  ^ 

ill 


=  (Btu/lb)  latent  heat  of  vaporization  of  saturated 
water  vapor  at  Ts<| 

=  The  fraction  of  surface  body  area  A  that  is  complete- 

* 

ly  wet  with  sweat  (moisture  loss  from  the  body  area 
is  imagined  to  occur  as  though  the  body  were  divided 
into  two  distinct  conditions — one  portion  entirely 
wet  and  the  remainder  dry).  Alternatively,  x  can  be 
envisioned  as  a  fraction  of  the  sweat  glands  that  are 
active. * 

=  (lb  water  vapor/lb  dry  air)  humidity  ratio  that 
'  corresponds  to  saturated  air  at  Tg. 

=  (lb  water  vapor/lb  dry  air)  logarithmic  mean  average 
of  air  outlet  and  inlet  humidity  ratios  and  is  cal¬ 
culated  by  the  following  equation  (see  Appendix  2): 


wa  =  ws  - 


wo  -  wi 

w  -  W. 

ln 

s  o 


where  wq  and  w^  are  the  humidity  ratios  of  air  at  exhaust  and  in¬ 
let  conditions,  respectively. 

Upon  substitution  of  Eqs.  (3)  and  r5)  in  Eq.  (2): 

HR  =  h  A  (T  -  T  )  +  h  h-  Ax  (we  -  w  )  (7) 

c  s  a  m  fg  s  a 

In  order  to  apply  Eq.  (8)  to  predict  equal  comfort  states 
for  the  human,  one  must  determine  values  of  heat  and  mass  trans¬ 
fer  coefficients  and  devise  a  method  of  estimating  x  and  the 
corresponding  value  of  T  under  different  environmental  condi- 
tions. 

Strictly  speaking,  water  loss  is  from  respired  air,  passive 
secretion  from  the  skin  and  active  sweat  glands.  For  simplicity, 
the  respired  and  passive  skin  losses  are  considered  as  an  equiva¬ 
lent  fraction  of  active  sweat  glands. 
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Natural  Motion  Limit 

In  Reference  1,  the  natural  motion  limit  concept  for  con¬ 
vective  heat  and  mass  transfer  coefficients  was  established. 
Attempts  to  correlate  these  coefficients  from  experimental  ev¬ 
idence  produced  by  the  calorimeter  on  the  basis  of  natural  or 
forced  convection  parameters  proved  unsuccessful.  It  was  con¬ 
cluded  that  the  natural  body  movements  of  the  subject  were  in¬ 
ducing  boundary  air  motions  which  were  overriding  the  weaker 
fluid  forces  produced  by  buoyant  or  mechanical  means  (5-15  cfm/ 
person).  Hence,  for  the  very  low  air  flow  rates,  limiting  values 
of  the  heat  and  mass  transfer  coefficients  can  be  assumed  as 
constants.  Obviously,  when  the  buoyant  or  mechanical  induced 
fluid  forces  create  heat  and  mass  transfer  coefficients  that 
exceed  these  limits,  the  constancy  of  the  values  will  no  longer 
pertain . 

Evaluation  of  Convective  Heat  Transfer  Coefficient 

The  natural  motion  limit  convective  heat  transfer  coeffi¬ 
cient  hc  is  determined  directly  from  the  calorimeter  experiments 
by  using  the  following  equation: 


The  value  of  C  is  found  by  measuring  air  flow  rate  and 
change  in  sensible  air  state  between  inlet  and  exit.  The  average 
value  of  hc  for  all  tests  is  found  to  be  equal  to  a  .76  Btu/hr 
ft2°F. 


Evaluation  of  Convective  Mass  Transfer  Coefficient  and  Per  Cent 


In  the  sweating  region  when  the  skin  is  not  completely  wet, 
it  is  very  difficult  to  compute  precise  values  for  h^  and  x. 


There  are  no  reliable  direct  methods  or  instruments  known  to  the 
authors  to  measure  the  value  of  x.  An  attempt  was  made  to  eval¬ 
uate  hm  by  asking  test  subjects  to  estimate  the  wetness  of  their 
bodies  from  which  a  determination  of  x  could  be  made.  Evalua¬ 
tion  of  many  test  results  [8]  revealed  that  this  method  was 
highly  inaccurate  as  was  evident  from  data  scatter  and  as  will  be 
discussed  later,  resulted  in  an  overestimation  of  x. 

An  improved  technique  is  involved  in  the  present  study. 
First,  the  h^x  product  is  determined  and  then  plotted  against 
average  skin  temperature,  Tg.  The  value  of  hmx  can  be  readily 
determined,  independent  of  test  subject  response,  with  reasonable 
accuracy  directly  from  the  calorimeter  data.  Since  hm  is  dic¬ 
tated  by  test  conditions  described  as  natural  motion  limit  condi¬ 
tions,  hm  is  assumed  constant  for  all  tests.  The  variation  of 
the  product  hmx  with  Tg  will  then  represent  only  the  variation  of 
x  until  the  whole  body  is  completely  wet  and  x  becomes  unity. 
Beyond  this  point,  x  remains  constant,  the  human  body  is  no 
longer  able  to  compensate  for  atmospheric  conditions  by  sweating 
and  the  deep  body  and  skin  temperature  must  rise  until  symptoms 
of  heat  stroke  appear. 

By  plotting  the  values  of  h  x  versus  T  for  all  tests  and 

ms 

fitting  curves  to  this  data  the  h  x  -  T  relation  shown  in 

ms 

Fig.  2  is  obtained.  The  asymptotic  horizontal  line  gives  the 
value  of  h  since  x  is  unity  here.  The  natural  motion  limit 

2 

value  of  h^  is  found  from  Fig.  2  to  be  1.4  Ibm  water  vapor/hr  ft 

/ 

(lbm  water  vapor/lb  dry  air).  With  the  value  of  hm  known,  the 

relation  of  x  only  versus  T  can  be  established,  as  shown  in 
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FIGURE  2 
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BOOY  WETNESS  FRACTION  VS.  SKIN  TEMPERATURE 


X<  FRACTION  OF  SKIN  SURFACE  THAT  IS  WET 
T*j  BODY  AREA  AVERAGED  SKIN  TEMPERATURE  #F 


85  86  87  88  89  90  91  92  93  94  95  96  97  98  99  100  101  102 


FIGURE  3 
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III.  THE  RELATIONSHIP  BETWEEN  (x  -  T  ) 

AND  STATES  OF  EQUAL  COMFORT 

Studies  have  been  made  of  the  behavior  of  evaporation  rates 
from  the  human  body  under  various  environmental  conditions.  One 
of  these  studies  was  by  C.  E.  Winslow,  0.  P.  Herrington  and 
A.  P.  Gagge  [113.  They  found  the  total  evaporative  heat  loss  is 
independent  of  atmospheric  humidity  but  dependent  on  air  tem¬ 
perature.  This  can  be  explained  by  physiological  control  of 
sweat  secretion.  If  the  surface  of  the  body  maintains  the  same 
physiological  characteristics,  i.e.,  maintains  the  surface  wet¬ 
ness  and  skin  temperature,  there  would  be  less  evaporation  with 
high  relative  humidity.  The  rate  of  evaporation  from  a  unit 
area  of  wet  surface  is  proportional  to  the  difference  between 
water  vapor  concentration  at  the  surface  and  that  of  the  sur¬ 
rounding  air  providing  air  movement  is  constant*  The  constancy 
of  evaporative  heat  loss  with  varying  relative  humidity  must  be 
maintained  by  change  in  sweat  secretion.  As  the  humidity  ratio 
of  the  surrounding  environment  is  increased,  there  occurs  within 
a  given  area  an  increase  in  the  number  of  active  sweat  glands , 
and  additional  glands  are  recruited  in  previously  nonsweating 
areas.  This  increase  in  the  number  of  the  active  sweat  glands 
continues  with  more  severe  environments  until  the  limit  is 
reached*,  the  recruitment  of  sweat  spreads  to  all  skin  surfaces. 

It  should  be  noted  from  Fig.  3  that  there  is  a  slight  rise 
in  skin  temperature  as  the  sweat  glands  are  recruited.  Since 
this  results  in  a  small  reduction  in  convective  heat  transfer, 
the  evaporative  heat  loss  will  increase  slightly  as  the  air  humid¬ 
ity  increases.  From  a  strict  point  of  view,  the  finding  of 
Winslow,  et  al. ,  is  not  correct. 


.n.„  .  . 
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Interpretation  of  (x  -  Tg)  Relationship 

In  Section  II  the  relation  between  skin  wetness,  x,  as  the 
fraction  of  the  body  surface  that  is  covered  by  an  evaporative 
thin  film  of  sweat,  and  the  average  skin  temperature  was  estab¬ 
lished;  see  Fig.  3.  This  is  the  principal  physical  property 
relationship  of  the  boundary  skin  surface  and  is  the  resultant 
of  both  various  internal  regulatory  processes  and  the  external 
stress  caused  by  the  environment.  Therefore,  every  point  on  the 
(x  -  T  )  curve  of  Fig.  3  is  a  unique  physiological  comfort  or 


discomfort  condition.  If  the  human  body  is  exposed  to  different 


environmental  conditions  that  would  keep  the  value  of  T_  and  the 


corresponding  value  of  x  constant  irrespective  of  the  heat  ex 


change  mechanism,  then  these  different  environmental  states  can 


be  defined  as  states  of  equal  comfort. 


Subdivision  of  the  (x  -  Tg)  Relationship 

The  (x  -  T  )  relation  can  be  divided  into  three  distinct 
s 

regions : 

( 1)  Insensible  or  Passive  Sweating  Region.  This  region  is 
bounded  by  a  skin  temperature  below  94° F  and  above  the  shivering 
point.  The  skin  wetness  is  almost  constant  due  to  the  absence  of 
sweat  secretion.  The  skin  temperature  is  the  dominant  thermal 
regulatory  control,  i.e.,  the  body  temperature  regulation  is 
mainly  by  change  in  skin  temperature.  There  is  some  uncertainty 
as  to  whether  or  not  there  is  an  exact  value  of  the  critical  skin 
temperature  at  which  sweat  secretion  starts.  This  is  due  to  the 
lack  of  sufficient  tests  around  this  critical  point.  Hence,  the 


transition  between  insensible  and  sensible  sweat  regions  is  not 
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known  precisely  as  to  whether  it  is  a  smooth  or  abrupt  transition. 
There  is  some  evidence  based  on  tests  of  one  individual  to  support 
the  view  that  the  transition  is  smooth)  i.e.,  in  the  transition 
region  active  sweating  is  conditioned  by  gradual  rise  in  skin 
temperature  after  which  the  rate  of  increase  of  sweat  secretion 
with  skin  temperature  raches  a  maximum  value. 

(2)  Modulated  Active  or  Sensible  Sweating  Region.  This  region 
is  bounded  by  Tg  between  94°F  and  96.2°F.  As  mentioned  above, 
in  this  region  the  rate  of  increase  of  wetness  with  skin  tem¬ 
perature  reaches  maximum  limit.  In  other  words,  the  rise  of  skin 
temperature  is  very  small  compared  to  the  increase  of  x.  This 
occurs  because  the  sweat  glands  suddenly  begin  to  secrete  and  are 
able  to  dissipate  heat  by  evaporation  so  successfully  that  no 
further  significant  rise  in  skin  temperature  occurs. 

(3)  Heat  Stroke  Region.  In  this  region,  Ts  >  96.2°F,  there  is 

a  loss  of  thermal  regulation.  The  body  surface  is  completely  wet 
and  the  body  is  no  longer  able  to  dissipate  more  heat  by  evapora¬ 
tion  even  though  the  sweat  secretion  increases.  Under  this  condi¬ 
tion  the  skin  temperature  starts  to  rise  in  a  fashion  similar  to 
region  1.  Since  no  further  body  surface  regulation  is  possible, 
the  deep  body  temperature  will  also  rise  until  the  body  exceeds 
high  temperature  tolerance  and  heat  stroke  begins. 

In  the  present  analysis,  only  regions  one  and  two  are  con¬ 
sidered,  since  it  is  in  these  two  regions  only  that  the  body 
maintains  normal  thermal  regulatory  equilibrium  with  the  environ¬ 
ment. 

In  summary,  the  principal  phenomena  observed  from  the 

physiological  relation  of  Cx  -  T_)  are: 

s 


f: 


k 


}rf 

¥ 

& 


L 

& 


I 

I 

I 

T 


18 


(a)  The  triggering  of  the  active  sweat  secreting  mechanisms  is 
intimately  related  to  a  critical  skin  temperature.  The 
average  value  of  this  critical  temperature  is  in  the 
vicinity  of  94° F. 

(b)  There  are  two  major  thermal  regulators  controls  on  the 
human  body  surface,  namely,  the  skin  temperature  and  sweat 
secretion. 

(c)  Points  on  the  x  -  Tg  graph  define  different  settings  of  the 
body's  thermoregulatory  function  and  hence  are  points  of 
thermal  comfort  or  discomfort. 

Physiological  Thermal  Index 

Values  can  be  assigned  to  points  on  the  (x  -  T^)  relation 
to  give  the  level  of  physiological  strain  or  level  of, external 
environmental  stress.  Any  value  from  the  physiological  scale,  a 
point  on  the  (x  -  Tg)  curve,  represents  a  fixed  physiological 
state  regardless  of  the  particular  environment  which  produces 
this  state.  Because  there  are  two  different  thermal  regulatory 
controls  on  the  human  body  surface,  one  the  sensible  heat  regula¬ 
tory  control,  Ts,  and  the  other  the  evaporative  heat  control,  x, 
a  reasonable  index  should  weigh  these  two  thermal  regulatory 
controls  equally  according  to  their  power  to  dissipate  the 
metabolic  heat  output  under  the  same  environmental  condition. 

The  method  proposed  to  evaluate  the  heat  dissipating  power  by 
varying  Tg  as  compared  to  varying  x  is  to  calculate  the  ratio  of 
metabolic  heat  output  dissipation  by  unit  change  in  T  and  assign 
a  unit  to  a  change  in  x  that  will  result  in  an  equal  heat  dissi¬ 
pating  power.  From  this  heat  dissipating  power  ratio,  increments 


#< i&A^gptf* «*M^jafc3fcs  «sh;  -  -g£^  -  v 


of  x  and  T_  can  be  weighed  accordingly  on  the  scale.  The  heat 

5 

dissipating  power  ratio  is  derived  as  follows: 

Under  constant  environmental  conditions,  the  metabolic 

rate  is  a  function  of  T  and  x.  This  function  can  be  written  as: 

s 


&ij  »*» 


I!  I 


MR  a  f  (Tg ,  x)  (9) 

This  relation  shows  tl  a c  if  metabolic  heat  output  of  a  human 
body  is  incrementally  shifted  while  holding  a  constant  environ¬ 
ment,  the  only  way  to  dissipate  this  increased  amount  of  heat 
output  is  by  change  in  the  values  of  T^,  and  x.  This  incremental 
change  in  metabolic  rate  can  be  expressed  as  follows: 


dMR  *  c«-)xaTs  ♦  (||)T  dx 

S  S 


if  dMR  =  0. 
Therefore , 


From  Eqn.  (7)  , 


dx  _ 


/  3f  x 


~yT, 


(10) 


(11) 


CW~\  *  hCA 


‘jxV  =  hm  hfg  A  (ws  -  V 

s 


Upon  substituting  Eqs.  (12)  and  (13)  into  Eqs.  (11), 


a  h-  (we  -  wj 
m  fg  s  a 


(14) 


Equation  (14)  shows  is  dependent  on  air  humidity,  i.e.,  the 

Q1s 

cooling  power  ratio  of  evaporation  to  convective  heat  transfer 
varies  with  the  air  humidity.  The  evaporation  rate  per  unit 
surface  area  is  a  maximum  when  air  humidity  is  zero  and  it  is 


minimum  when  the  air  humidity  ratio  is  equal  to  that  at  the  skin 
surface.  The  proposed  index  is  built  by  choosing  a  value  for 
which  is  evaluated  at  wa  equal  to  zero.  As  a  first  approxi¬ 
mation,  the  right-hand  side  can  then  be  assumed  a  constant  with 
h^g  and  wg  evaluated  for  Tg  =  95°F  and  substituting  the  experi¬ 
mentally  determined  values  of  hc  and  h^  in  Eqn.  (14),  the  result 


art =  —  Tir (or)'1  (1S 

O 

Equation  (15)  shows  that  when  the  surrounding  environment  is 
completely  dry,  a  change  of  x  from  zero  to  completely  wet  body 
surface  is  equivalent  to  a  rise  of  70°F  in  skin  temperature. 

With  the  heat  dissipating  power  ratio  established,  the 
physiological  index  can  be  constructed.  The  zero  point  on  the 
index  is  taken  to  be  the  point  on  (x  -  Tg)  curve  of  skin  tem¬ 
perature  of  93. 4° F.  The  corresponding  body  temperature  under 
equilibrium  condition  is  found  to  be  97.9°F  [15].  The  values 
93.4°F  and  97.9°F  have  been  observed  as  the  average  temperatures 
of  the  skin  and  deep  body,  where  there  is  minimal  regulatory 
effort  in  maintaining  body  temperature  either  by  any  vascular 
effort  or  by  sweating.  When  these  temperatures  occur  simulta¬ 
neously  during  rest,  the  body  is  in  a  state  of  "Physiological 
Thermal  Nuetrality"  [14].  Above  this  natural  point,  the  body 
will  sense  warmth  and  the  reading  on  the  index  will  be  positive 
and  below  it  the  body  will  sense  coolness  and  the  reading  on 
the  index  will  be  negative.  The  index  is  divided  in  such  a  way 
that  one  division  represents  a  change  of  .05  in  x  or  3.5°F 
change  in  T_  or  combinatorial  change  of  T_  and  x  weighted  accord 
ing  to  their  heat  dissipating  power.  The  index  is  superimposed 


■  |  T 
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on  (x  -  T  )  relation  as  shown  in  Fig.  4.  Each  division  on  the 


index  represents  equal  amounts  of  metabolic  heat  output  that  can 


be  dissipated  by  physiological  change  of  x,  T  ,  or  both.  This 

s 


scale  is  defined  as  the  Physiological  Thermal  Index  (PTI). 


In  Section  IV  the  PTI  is  combined  with  the  environmental 


conditions  through  heat  and  mass  transfer  rate  equations  to 


predict  environmental  states  of  equal  comfort.  These  states  of 


equal  comfort  are  defined  as  the  air  state  that  would  result  in 


the  same  values  of  x  and  T  . 

s 


Numerous  problems  have  been  identified  and  numerous 


attempts  have  been  made  to  improve  the  effective  temperature  con¬ 


cept-,  some  of  which  have  been  cited  previously.  One  of  the  most 


serious  difficulties  associated  with  the  effective  temperature 


and  other  indices  of  comfort  is  that  they  are  tied  to  environmen¬ 


tal  conditions  as  part  of  their  basic  definition.  With  the 


Physiological  Thermal  Index  established  independently  of  envi¬ 


ronmental  conditions ,  it  becomes  the  primary  index  of  human 


comfort.  The  same  PTI  line  will  be  located  at  different  positions 


on  the  psychrometric  chart  depending  upon  metabolism,  air  flow 


rate,  temperature,  humidity  and  radiation  factors.  Yet,  the  line 


will  represent  the  same  index  value  on  the  x  -  Ts  relationship 


and  hence  a  constant  index  of  comfort. 
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PHYSIOLOGICAL  THERMAL  INDEX  (PTI)  SUPERIMPOSED  ON  X  VS.  Ts 
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IV.  THE  LOCI  OF  PTI  VALUES  ON  THE  PSYCHROMETRIC  CHART 

In  locating  lines  of  constant  PTI  values  on  the  psychro- 
metric  chart  the  following  simplifications  and  assumptions  are 
made : 

(1)  The  specific  heat  and  molecular  weight  of  dry  air  are  used 
in  place  of  the  water  vapor-dry  air  mixture  values,  since  the 
mixture  contains  no  more  than  3  per  cent  water  vapor  by  weight. 
(‘2)  The  total  evaporative  heat  losses  which  include  heat  of 
vaporized  moisture  from  the  lungs  during  respiration,  heat  of 
vaporized  water  diffusing  through  skin  layer  and  the  heat  of 
vaporized  sweat  necessary  for  the  regulation  of  body  tempera¬ 
ture  are  taken  as  one  quantity  and  represented  by  a  fraction,  x, 

,  of  the  skin  that  is  assumed  saturated. 

(3)  The  human  body  is  in  steady  state  with  the  environment, 
i.e.,  the  average  skin  temperature,  deep  body  temperature  and 
metabolic  rate  are  not  varying  with  time. 

(4)  The  humidity  ratio  at  the  wet  skin  surface  is  determined 
assuming  a  saturated  condition  at  skin  temperature,  T  . 

b 

1  Derivation  of  Psychrometric  Equations  Using  PTI  Values 

The  generated  lines  of  equal  comfort  are  based  on  average 
environmental  conditions.  In  the  human  calorimeter,  which  is 
a  relatively  confined  space,  or  in  shelters  with  minimum  ventila¬ 
tion,  the  test  subjects  (due  to  their  energy  and  moisture 
release)  will  experience  significant  differences  between  the 
1  initial  air  state  as  it  contacts  their  body  and  final  air  state 
a,s  it  leaves  their  body.  For  this  reason  the  logarithmic  average 
of  inlet  and  exhaust  states  is  assumed  as  the  average  air  condi- 
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tion  to  which  the  test  subject  is  exposed.  Sensible  and  evapora¬ 
tive  heat  outputs  and  average  air  states  are  calculated  in  terms 
of  the  inlet  and  exit  air  states  by  the  following  equations: 

Sensible  heat  is  calculated  in  terms  of  air  mass  flow  rate, 
inlet  and  exhaust  dry  bulb  temperature  by 

c  =  G  cp  <To  "  Ti)  (16) 

where  G  (lbm/hr)  is  the  air  mass  flow  rate.  C  (Btu/lbm°F)  is 

P 

the  specific  heat  of  dry  air  at  constant  pressure. 


From  Eqs.  (3)  and  (4),  C  can  be  rewritten  as: 


C  =  h„  A 
c 


In  ( 


."s 

T” 
s  o 


Combining  Eqs.  (16)  and  (17),  Tq  can  be  written  in  terms  of  as 


follows : 


where 


To  =  Ts  (1  -  e'Y)  +  T.  e~Y 


.  hc  A 

‘  '  s-c; 


(18) 


Substituting  Eq.  (18)  into  Eq.  (17),  C  can  be  recast  in  terms  of 
inlet  dry  bulb  temperature,  mass  flow  rate  and  convective  heat 


transfer  coefficient,  as: 


C  =  G  Cp  (Tg  -  T£)  (1  -  e"1) 


(19) 


The  total  evaporative  heat  loss  can  also  be  related  to 
humidity  ratios  at  inlet  and  exit  states  and  air  mass  flow  rate 
by  the  following  equation: 


E  =  G  hfg  (wQ  -  w^) 


(20) 
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Upon  substitution  of  Eq.  (6)  into  Eq.  (5),  E  can  be  rewritten  as: 


Combining  Eqs. 
as  follows: 


E  =  h  h-  Ax 
m  fg 


wo  -  w. 


w  -  w. 
In  (~ - ~) 


(21) 


w  -  wrt 
s  o 


(20)  and  (21),  w^  can  be  written  in  terms  of  w. 

o  1 


wrt  =  w„  -  (w„  -  w. ) 
os  si 


-Z 


(22) 


where 


Z 


h  Ax 
m 

g 


Substituting  Eq.  (22)  into  Eq.  (20) 

E  =  G  hfg  (ws  -  w£)(l  -  e"Z)  (23) 

From  Eqs.  (23),  (22)  and  (6),  total  evaporative  heat  loss,  humid¬ 
ity  ratio  at  exit  condition  and  average  humidity  ratio  can  be 
determined  for  constant  skin  temperature  and  constant  wetness , 
but  at  varying  air  mass  flow  rate,  by  substituting  different 
values  of  w^  in  the  corresponding  equations. 

Sensible  heat  output,  C,  can  be  determined  by  subtracting 
the  evaporative  heat  loss  from  the  metabolic  heat  output.  Know¬ 
ing  the  sensible  heat  output, mass  flow  rate,  and  skin  temperature, 
the  inlet  dry  bulb  temperature  can  be  calculated  by 


T.  =  Ts  - 


G  Cp  (1  -  e”Y) 


(24) 


Note:  The  preceding  relations  are  applied  to  calculate  the  inlet 
and  exit  air  states  under  various  air  flow  rates.  The  average 
air  state,  sensible  and  evaporative  heat  output  can  be  held 
constant  while  varying  inlet  states  and  the  air  flow  rates. 
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A  computer  program  is  written  (see  Appendix  3)  to  compute 
average  air  states,  sensible  and  evaporative  heat  outputs,  inlet 
and  exit  air  states  that  correspond  to  constant  skin  temperature 
and  constant  wetness,  i.e.,  air  states  of  equal  comfort.  Meta¬ 
bolic  rate  may  also  be  varied  as  a  program  parameter. 

Generation  of  Lines  for  PTI  Constant  Values 

Two  sets  of  PTI  constant  lines  are  generated  by  the  comput¬ 
er  program.  One  set,  shown  in  Fig.  5,  represents  a  typical 
population  sample  (based  on  1960  United  States  Census  [12])  which 

has  an  average  metabolic  rate  of  275  Btu/hr  Cl]  and  surface  area 
2 

of  15.6  ft  .  The  other  set,  shown  in  Fig.  6,  represents  adult 
males  with  average  metabolic  rate  of  400  Btu/hr  [13]  and  surface 
area  of  20  ft2.  C.  P.  Yaglou’s  ET  lines  for  nude  and  clothed 
subjects  [17]  are  drawn  in  both  Figs.  5  and  6,  as  comparisons  to 
the  predicted  PTI  constant  lines. 

Characteristics  of  Predicted  PTI  Lines 

The  predicted  PTI  lines  (loci  of  constant  skin  temperature 
and  constant  wetness)  are  shown  as  solid  lines  in  Fig.  5.  The 
constructed  PTI  lines  cover  the  range  of  slightly  cool  environ¬ 
ment  (-3.4  PTI  line)  to  the  maximum  tolerable  limit  where  the 
body  is  completely  wet.  This  maximum  limit  is  represented  by 
the  17.64  PTI  line.  It  is  noticed,  in  the  slightly  cool  environ¬ 
ment  that  the  PTI  lines  are  very  steep.  As  the  PTI  value  in¬ 
creases  the  slope  of  the  line  decreases  to  a  minimum  value  which 
is  represented  by  the  17.64  PTI  line.  Beyond  this  PTI  line  the 
slope  remains  the  same  because  the  surface  wetness  remains 
constant  (x  =  1)  and  the  skin  temperature  rises. 
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From  this  observation,  it  is  concluded  that  the  humidity 
ratio  of  the  environment  has  slight  effect  on  comfort  at  lower 
PTI  values  while  in  the  higher  PTI  range  it  has  a  pronounced 
effect.  This  can  be  interpreted  by  noting  that  the  skin  tem¬ 
perature  is  the  dominant  index  of  comfort  in  the  region  of  in¬ 
sensible  sweating  while  the  skin  wetness  is  the  dominant  index 
of  comfort  or  discomfort  in  the  active  sweating  region. 

C.  P.  Yaglou*s  ET  lines  are  drawn  in  Fig.  5  in  broken 
lines  as  a  comparison  to  the  predicted  PTI  lines.  But  we  cannot 
draw  a  comparative  conclusion  between  these  lines  since  Yaglou* s 
ET  lines  are  for  adult  males  with  different  metabolic  rates  per 
unit  of  body  area  than  the  general  population.  Although  the 
generated  PTI  lines  for  adult  males  have  the  same  general  trends 
as  that  of  the  general  population  except  as  will  be  discussed 
later  it  is  more  appropriate  to  limit  the  comparison  only  between 
ET  lines  for  nude  adult  males  generated  by  Yaglou  and  by  our 
predictive  model. 

Comparison  Between  Yaglou1 s  Values  and  Predicted  PTI  Lines  for 
Adult  Males 

A  comparison  between  Yaglou  ET  lines  and  PTI  lines  pre¬ 
dicted  by  the  present  study  can  be  made  by  superimposing  each 
set  of  ET  lines  on  the  same  psychrometric  chart,  as  shown  in 
Fig.  6.  The  predicted  PTI  lines  are  shown  by  solid  lines  while 
Yaglou  lines  are  shown  by  broken  lines.  In  the  region  of  lower 
PTI  lines,  it  is  noticed  that  there  is  a  great  deviation  between 
Yaglou  lines  and  the  predicted  ones.  Yaglou  lines  show  a  much 
greater  humidity  dependence  than  the  predicted  lines.  In  the 
region  of  higher  PTI  lines,  Yaglou  lines  correspond  quite  well 
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with  the  predicted  lines,  particularly  at  the  85  ET  line.  On 
the  other  hand,  the  predicted  line  of  x  equals  unity  shows 
greater  dependence  on  humidity  than  the  corresponding  Yaglou  ET 
line.  From  the  preceding  analysis,  we  can  draw  the  same  conclu¬ 
sion  as  other  workers  that  Yaglou  ET  lines  overestimate  the  in¬ 
fluence  of  humidity  on  sensations  of  warmth  and  comfort  at  or¬ 
dinary  temperatures  and  underestimate  the  effect  at  very  high 
temperatures. 

This  conclusion  can  also  be  arrived  at  by  examining  the 
variation  of  wetness  (x)  with  average  dry  bulb  temperature  using 
the  Yaglou  ET  line  as  a  parameter— see  Fig.  7.  The  wetness  x  at 
any  dry  bulb  temperature  along  the  Yaglou  ET  line  is  calculated 
by  the  following  equation: 


x 


MR  -  h  A  (T6  -  TJ 
_ c  s  a 

h  h  'A  (w  -  w  ) 
m  rg  s  a 


(25) 


where 

MR  is  the  average  metabolic  rate  output  of  adult  males 
(400  Btu/hr) . 

hc  is  the  experimentally  determined  convective  heat 

transfer  coefficient  (.76  Btu/hr). 

A  body  surface  area  for  average  adult  male  (20.0  sq  ft). 

hm  Experimentally  determined  convective  mass  transfer 

2 

coefficient  (1.4  lb  water  vapor/hr  ft  (lb  water 
vapor/ lb  dry  air)). 

hfg  latent  heat  of  vaporization  at  the  partial  pressure  of 
saturated  water  vapor  at  T  . 

Ta  air  dry  bulb  temperature  at  a  point  on  Yaglou  ET  line 
that  corresponds  to  T  . 

a 


X’  FRACTION  OF  BOOY  SKIN  THAT  IS  WET 

Te*  BODY  AREA  AVERAGED  SKIN  TEMPERATURE  °F 
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THERMAL  COMFORT  ZONES  ON  (X-T)  CURVE 


FIGURE  8 


COMFORT  ZONES  FOR  SEMI-NUDE  ADULT  MALES 
UNDER  VERY  LOW  AIR  FLOW  RATES 
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Fig.  10.  Our  predictive  PTI  lines  and  F.  Rohles  and  R.  Mevins 
comfort  lines  [22]  are  also  drawn  in  Fig.  10  for  comparison.  It' 
should  be  pointed  out  that  Nevins  and  Gagge  used  different  valur' 
environmental  parameters  than  we  used.  The  present  study  deals 
primarily  with  semi-nude  people  under  the  environment  of  very 
low  ventilation  rates  and  in  the  absence  of  radiation  exchange , 
while  Gagge  and  Nevins  used  the  environmental  characteristics  of 
lightly  clothed  people  and  wall  temperatures  equal  to  air  tem¬ 
perature.  Therefore,  the  predicted  PTI  lines  and  Gagge’ s  new  ET 
lines  that  coincide  do  not  necessarily  represent  the  same  Uv.n.. 
of  comfort.  But  the  important  observation  noticed  from  Fig.  1G 
is  that  the  PTI  lines  have  the  same  slopes  as  those  of  Gagge,  i.~-.. 
they  agree  well  on  humidity  dependence  and  on  thermal  sonsatir  n . 
Subjective  Experimental  Response  Compared  to  the  Hew  PTI  Lines 
At  the  time  the  new  PTI  lines  were  being  developed  some 
tests  in  the  lower  temperature  range  were  performed  on  subjects 
in  the  monoman  calorimeter  with  each  test  consisting  of  two 
parts.  In  the  first  part,  the  inlet  dry  and  wet  bulb  temperatures 
were  chosen  and  kept  constant  during  that  part  of  the  tes.  and 
data  was  collected  after  the  subject  reached  equilibrium  with 
his  environment.  The  second  part  of  the  test  consisted  of 
keeping  the  dry  bulb  temperature  the  same  as  that  of  the  first 
part  of  the  test  and  changing  the  wet  bulb  temperature  to  a 
different  value.  After  the  subject  reached  equilibrium  in  the 
second  part  of  the  test,  verbal  expressions  were  obtained  from 
the  subjects  regarding  their  sense  of  hotness,  coolness  or  com¬ 
fort  to  the  first  state.  Surprisingly,  most  of  the  answu  r>  om 
the  test  subjects  were  that  they  sensed  the  same  comfort  as  that 


DRY  BULB  TEMPERATURE  °F 


of  the  first  part  of  the  test  even  though  the  majority  did 

observe  the  change  in  humidity.  From  these  tests  it  appears'  , 

that  the  humidity  has  little  effect  on  comfort  in  the  normal  ten  * 

perature  range.  This  supports  the  new  method  for  prediction  of 

;  ■ 

the  PTI  lines.  Although  the  number  of  these  tests  is  limited  to 
date,  they  give  significant  support  to  the. concepts  underlying 
the  physiological  thermal  index.  '  '  . 

Work  conducted  during  this, project  but  previous  to,  the 

} 

creation  of  the  predictive  model  involving  the  PTI  index  lends 
further  support  to  the  worth  of  the  concept.  .  1 

Comparisons  of  psychological  impressions  with  physical  . 

■  i 

environment  were  attempted  in  32  pairs  of  tests  in  which  th<= 
test  subject  was  exposed  sequentially  to’ two  different  climate 

a  *  » 

states.  The  test  conditions  were  82  and  85°F  ET  ait  inlet  state 

s 

ana  13  cfm.  Subject  thermal  equilibrium'  was  established  for  each 

i 

test  state.  These  sequential  states  were  usually  on  the,  same 

*  : 

effective  temperature  line.  There  we're  21  pairs  of  tests  wit}. 

1 

constant  effective  temperatures  anid  11  pairs  of  tests  with 

,  I 

differing  effective  temperatures.  The  high  and  low  humidity 

sequences  were  alternated  from  one'  pair  of  tests  to  the  next. 

After  experiencing  , the  two  climate  states.,  the  subject  was 
asked  to  check  one  item  in  the  following  statement:’ 

"How  does  the  new  condition  compare  with  the’  first  condi¬ 


tion  of  about  ah  hour  ago?  More  comfortable 


Lefes 


comfortable 


Same 


In  effect,  the  subject  would  have  three  equally’ possible 
responses  if  the  response  were  random  compared  to  a  measured 
physical  parameter.  That  is,  the  most  probable  number  of 
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responses  of  one  classification  to  the  total  number  of  responses 
would  be  one- third  if  there  were  no  correlation. 

A,  Response  Compared  to  Test  Sequence 

Those  who  indicated  1st  state  was  less  comfortable  7/32  2?% 
Those  who  indicated  2nd  state  was  less  comfortable  8/32  25% 


Those  who  indicated  the  states  were  the  same  17/32  5  "■ 

B.  Response  Compared  to  ET 

Those  whose  impression  agreed  with  the  change  in  21/32  ff° 
effective  temperature 

C.  Response  Compared  to  ql  (latent  heat)/qt  (total  heat) 


Those  who  felt  warmer  for  higher  ql/qt  5/32  16% 

Those  who  felt  cooler  for  higher  ql/qt  10 / ?2  3V. 

Those  who  felt  the  same  for  higher  ql/qt  17/32  53% 

D.  Response  Compared  to  Skin  Temperature* 

Those  who  were  less  comfortable  w5.th  higher  skin  7/2  8  2  5% 

temperature 

Those  who  were  more  comfortable  with  higher  skin  6/28  21% 

temperature 

Those  who  felt  the  same  with  higher  skin  15/28  25% 

temperature 
Plausible  Conclusions 

1.  Test  sequence  is  not  a  significant  factor  in  the  test  subject's 
response. 

2.  Effective  temperature  is  strongly  related  to  psychological 
impressions  for  these  test  conditions. 


17/32  5  3% 


7/2  8  2  5% 


6/28  21% 


15/2  8  2  5% 


The  reason  for  28  tests  here  is  that  some  subjects  produced 
the  same  skin  temperature  +  0.1°F  which  is  within  the  uncertainty 
estimate  of  the  average  skin  temperature. 
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3*  The  psychological  response  shows  no  significant  correlation 
with  changes  in  humidity  even  though  the  value  of  ql/qt  was 
changed  in  all  but  one  test  and  as  much  as  30%  in  some  ■■ases, 
Likewise,  no  significant  correlation  is  observed  between 


K 


II 
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change  in  averaged  skin  temperatures  and  psychological 
response.  It  is  to  be  noted  that  skin  temperature  changes 
for  the  two  test  states  were  usually  less  than  1/2°F. 

It  was  shown  in  Fig.  7  for  the  82-B5°F  effective  tempera¬ 
ture  range  that  there  was  good  agreement  between  the  PTI  lines 
and  the  ET  lines.  The  experimental  evidence  just  cited  gives 
physical  testimony  as  to  the  correctness  of  the  predictive  model 
in  this  region  and  further  strengthens  the  merits  of  the  Physio¬ 
logical  Thermal  Index  concept. 

Predicted  PTI  Lines  of  Adult  Males  Compared  with  General 
Population 

The  PTI  lines  for  a  general  population  sample  and  those 
for  adult  males  are  superimposed  in  Fig.  11.  It  is  in:eresting 
to  note  in  this  figure  that  the  slopes  of  PTI  lines  that  lie 
on  the  psychrometric  chart  in  the  range  of  65-70°F  saturation 
temperature  for  both  adult  males  and  the  general  population  are 
nearly  the  same  but  with  significantly  different  values  of  x  and 
T  .  This  means  that  the  effect  of  humidity  on  comfort  for  both 

O 

adult  males  and  the  general  population  is  the  same  in  this  range 

of  PTI  lines.  The  different  values  of  T  for  adult  males  and 

s 

general  population  at  65° F  saturation  indicates  that  they  should 
not  sense  the  same  comfort  under  the  same  environment  due  to  the 
difference  i<.  their  metabolic  rate  per  unit  of  body  surface  area. 

At  higher  PTI  lines  adult  males  show  a  greater  dependence  on  humid- 
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ity  for  their  comfort  than  the  general  population.  Assuming 
equal  acclimatization  and  tolerance  the  general  population  should 
be  able  to  tolerate  more  severe  environmental  conditions  than 
adult  males,  as  is  evidenced  by  contrasting  the  PTI  lines  at  >: 
equals  unity  for  adult  males  and  the  general  population. 

The  preceding  discussion  indicates  that  subjects  with 
different  metabolic  rates  and  different  body  areas  may  not  have 
the  same  physiological  comfort,  as  their  PTI  values  will  be  dif¬ 
ferent  under  the  same  environmental  conditions. 

Predicted  PTI  Lines  for  Subjects  in  Environments  of  High  Air 


It  is  concluded  from  the  previous  analysis  for  nude  seden¬ 
tary  subjects  at  thermal  equilibrium  for  very  low  air  ventilation 
rates  that  the  sense  of  equal  comfort  is  governed  by  the  constancy 
of  skin  temperature  and  wetness.  This  basic  principle  of  human 
comfort  can  be  extended  to  an  environment  with  different  air 
velocities  and  with  radiant  heat  exchange.  To  show  the  effect 
of  increased  air  velocity  and  radiant  exchange  between  the  human 
body  and  the  surroundings,  one  PTI  line  was  placed  in  three 
locations  on  the  psychrometric  chart,  as  shown  in  Fig.  12. 

The  solid  line  represents  the  comfort  state  for  the  general 
population  (T  =  85,  x  =  .134,  or  PTI  -  -  3.5)  for  very  low  air 

O 

ventilation  rate  and  no  radiation  exchange.  The  broken  line 
represents  the  same  sense  of  comfort  (same  x  and  Tg)  as  the 
solid  but  with  radiant  heat  loss  from  the  body  surface  to  the 
surroundings.  The  dotted  broken  line  represents  the  same  sense 
of  comfort  as  the  other  two  lines  but  with  high  air  velocity  and 
no  radiation  effects.  Because  the  heat  and  mass  transfer  coef- 
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ficients  are  dependent  upon  the  same  parameters,  by  increasiiIg 
the  air  velocity  and  the  heat  and  mass  transfer  coefficients  will 
increase  but  the  ratio  between  them  will  remain  constant.  It  is 
noticed  from  Fig.  12  that  by  increasing  the  air  velocity  the 
same  comfort  state  must  move  to  a  higher  environmental  condition 
but  the  line  of  comfort  of  high  air  velocity  is  parallel  to  the 
one  with  low  air  velocity.  With  the  presence  of  radiation  loss 
the  comfort  line  moves  to  a  higher  environment  but  with  steeper 
slope  than  that  without  radiation.  This  is  because  the  ratio  of 
sensible  heat  to  metabolic  heat  output  is  changed  with  the  pres¬ 
ence  of  radiation.  From  the  present  analysis  it  is  evident  that 
the  human  body  at  any  state  of  comfort  or  discomfort  can  be 
represented  by  a  single  point  on  the  (x  -  T  )  relation  while 
this  same  comfort  state  will  be  placed  differently  on  the  psy- 
chrometric  chart  depending  on  the  heat  and  evaporative  transfer 
mechanisms  and  air  velocity. 

Non-isostate  Process  Lines 

The  environmental  conditions  under  study  are  for  very  low 
air  flow  rates  and  confined  space.  The  air  flow  rates  range 
from  approximately  5  to  20  cfm,  which  give  negligible  mechanical¬ 
ly  induced  air  velocities  in  the  range  of  1  to  2  ft /min.  The 
average  FTI  line  will  vary  with  variation  in  the  air  flow  rate 
while  keeping  the  inlet  conditions  constant.  To  maintain  the 
same  average  PTI  line  the  air  state  at  the  inlet  condition  should 
be  changed  with  change  in  air  flow  rate.  Figures  13  and  14  show 
the  non-isostate  process  lines  for  the  average  PTI  lines  equiva¬ 
lent  to  ET  values  of  70  and  82. 5° F,  respectively.  The  non-isostate 
line  represents  loci  of  inlet  and  exit  states  for  different  air 
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flow  rates  that  would  give  the  same  average  physiological  thermal 
index.  It  is  noted  in  Fig.  13  that  the  exit  states  asymptotical¬ 
ly  approach  the  95° F  dry  bulb  temperature  which  represents  the 
skin  temperature  at  the  average  PTI  value  of  (82.5°F  ET). 
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V.  RESPONSE  TO  A  TRANSIENT  THERMAL  FLUX 


Review  of  Literature  Predictive  Models 

One  of  the  objectives  of  the  present  study  is  the  develop¬ 
ment  of  a  predictive  thermal  model  of  a  human  being  in  a  shelter 
environment  upon  exposure  to  transient  thermal  radiation  which 
simulates  the  radiation  emitted  bv  a  ceiling  with  fire  above  it. 

Two  models  from  the  literature  that  are  significant 
relative  to  this  program  objective  are  those  of  Machle  and  Hatch 
[23]  and  Kerslake  and  Waddell  [24]. 

The  Machle  and  Hatch  model  was  developed  beginning  with  a 
basic  energy  balance  equation  similar  to  Eq.  (1)  of  this  report 
(see  Section  II) 

MR-S=E  +  R+  C  +  W 
This  equation  was  then  expanded  as  follows: 

MR’  -  S  a  K  A  (T  -  T)  +  K  A  (T  -  T  )  +  K  AX  CVPc  -  VP) 
rsw  c  s  a  e  s  a 

where 

MR’  =  MR  -  W  (metabolic  rate  -  work  output) 

A  =  body  surface  area 
T  =  average  skin  temperature 

S 

Tw  =  mean  radiant  temperature  of  surroundings 
T_  s  air  temperature 

VP  =  vapor  pressure  of  water  at  temperature  T 

5  5 

VP  =  vaDor  pressure  of  moisture  in  the  air 

a 


(=  %  RH  x  VP  @  T  ) 

a 

X  =  evaporative  area/body  area 
Kr>K,and  Ke=  coefficients  of  heat  exchange  bv  radiation, 
convection  and  evaporation 
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The  storage  term,  S,  was  then  expanded  to 


dT. 


dT 


S  =  cm  (a 


dx 


r  +  bar> 


where 


c  = 
m  = 


dT  dT 
r  s 


dt  *dT 


average  body  specific  heat 
body  weight 

rates  of  change  of  average  values  for  rectal 
temperature  and  skin  temperature,  respectively. 
a,b  =  weighing  factors  to  proportion  the  body  heat 
storage  relative  to  the  deep  body  and  skin 
tissues 


Also,  Machle  and  Hatch  interlocked  T  and  T  with  the 

s  r 


relationship 


Tr  =  no  *  nl  Ts 


or 


dTr  =  n-L  dTs 


where  tj  ,  n-^  are  constants  determined  from  experimental  data. 

They  were  then  able  to  solve  for  the  percentage  of  equilib¬ 
rium  versus  time.  The  resulting  equation  is 


T  -  T 
eq  s 

TV  -"T 


exp  [- 


AC’t 


eq 


cimTTan1""+"'b7 


•3 


where 


C  :  K  +  K  +  2.5  K 


Teq=  new  equilibrium  skin  temperature 


TQ  =  initial  skin  temperature 


In  this  solution  it  was  assumed  that  M,  K  ,  T  ,  K  ,  T  ,  K  ,  X, 

X  W  U  Q  6 


and  VP=  were  constant  for  any  given  situation. 

a 


nil* 

-  ,  - 


m 


Upon  examination  of  this  equation  it  can  be  seen  that  the 
time  required  to  reach  a  given  percentage  of  the  new  equilibrium 
is  independent  of  both  the  surrounding  air  temperature  and  the 
humidity.  As  will  be  shown  later  this  is  not  believed  to  be 
factual;  the  response  time  is  dependent  upon  surrounding  air 
temperature  and  humidity. 

Another  problem  with  this  model  is  that  the  rectal  and 
skin  temperatures  are  locked  into  a  direct  functional  relation¬ 
ship  independent  of  any  other  parameters.  This  can  be  seen  as  a 
basic  weakness  except  for  very  slow  changes  in  environment.  This 
is  so,  for  the  skin  tissues  have  only  a  small  thermal  capacitance 
hence  they  can  respond  abruptly  to  changes  in  environmentally 
produced  heat  fluxes-  and  may  even  go  higher  than  the  rectal 
temperature  for  a  short  period  of  the  transient  response  time. 

As  explained  later  the  predictive  model  to  be  proposed  gets 
around  this  difficulty. 

Although  there  are  some  inadequacies  in  the  Machle  and  Hatch 
model  it  still  has  some  good  points.  It  appears  to  be  a  satis¬ 
factory  model  in  the  region  where  there  is  no  active  sweating, 
as  well  as  in  the  region  where  the  skin  is  completely  wot.  The 
model  might  have  been  expanded  to  include  the  active  sweat  region 
with  some  follow  on  work. 

Machle  and  Hatch  realized  in  their  study  that  an  environment 
which  would  produce  the  same  skin  temperature  as  a  second  environ¬ 
ment  could  be  considered  equivalent  to  the  second*  environment. 

They  were  handicapped,  however,  by  not  knowing  the  parameters 
which  influence  the  variation  of  skin  wetness. 
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Kerslake  and  Waddell  used  the  same  model  as  Machle  and  Hatch 
with  slightly  different  constants.  They  investigated  the  validitv 
of  the  model  when  the  skin  was  completelv  wet.  Thev  came  to  the 
conclusion  that  the  Machle  and  Hatch  equations  were  supported  bv 
their  experimental  results. 

Discussion  of  the  Transient  Thermal  Flux 

In  the  present  work,  a  predictive  technique  is  desired  for 
a  radiant  thermal  flux  which  simulates  a  building  fire  on  the 
floor  above  the  ceiling  of  a  survival  shelter.  The  literature 
[25]  provides  evidence  of  the  time  variation  and  intensity  of 
such  radiation  fluxes.  Several  buildings  were  burned  after 
being  instrumented  from  which  contrasting  Dulses  were  obtained. 

The  radiation  mode  chosen  for  the  present  study  was  that 
of  a  library  which  was  shown  to  have  the  highest  magnitude  of 
radiation  flux  of  any  of  the  buildings  burned.  The  air  tem¬ 
perature  and  humidity  were  assumed  to  remain  constant  during  the 
radiation  pulse.  Actually,  the  air  temperature  and  humiditv 
would  probably  not  stav  constant.  However,  to  obtain  a  reasonable 
basis  for  comparison  between  the  predictive  model  «-nd  experimental 
evidence  this  is  assumed  so.  If  there  is  good  agreement  between 
prediction  and  experiment,  no  great  difficulty  is  anticipated 
in  dealing  subsequently  with  transient  air  temperatures  and 
humidities . 


The  problem  then  is  to  predict  the  transient  thermal  behav¬ 
ior  of  a  human  subiect  in  terms  of  skin  and  deep  body  temperature 
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to  a  thermal  radiation  flux  characterized  by  a  library  fire  on 
the  floor  above  the  survival  shelter.  *••*  1  1 

1  1 

■  »  , 

The  assumptions  are: 

1 

(a)  The  sedentary  subject  is  in  thermal  equilibrium  with  his  1 

environment  prior  to  the  pulse.  ■  (  . 

(b)  He  maintains  his  metabolic  output  at  a  constant  level.  1  ■ 

1  •  J 

(c)  The  surrounding  air  temDerature  and  humiditv  remain  constant 

(d)  The  relationship  of  (x  -T  )  .developed  for  equilibrium 

b 

conditions  certains.  This  appears,  justified  by  the  fact 

«  1 

1  1  1 

that  the  response  of  the  sweat  glands  to  change  in  skin 

temperature  is  very1 fast  compared  to  the  response' of  the 

1 

core  and  skin  temperatures  to  the  proposed  change  in  .the 

1 

thermal  environment. 

5  t  *  1  S 

(e)  The  apparent  thermal  conductivity  between  the  skin  and  the 

« 

body  core  is  related  to  the  skin  temperature  as  reported 

’  1 

by  Winslow,  et  al.  [21].  In  their  studv,  the  reciprocal  , 

1  .1 

to  the  conductivity*  the  conductive  resistance  to  heat  flow, 
was  used.  Using  their  data,  Figw  15  was  constructed. 
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The  thermal  circuit  representing  the  subject  exposed  to 
the  proposed  transient  thermal  environment  is  shown  in  Fig.  16. 
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FIGURE  16 


MR  =  (Btu/hr)  metabolic  rate 


T  =  ( °F)  average  skin  temperature 
s 


(°F)  average  body  core  temperature — taken  as  a 


a,r 


weighted  mean  of  rectal  and  skin  temprature  [7] 

o  1 

T  =  £  T  +  £  T 
c  3  r  3  s 

=  (°F)  temperature  of  air  and  radiation  source, 


respectively 


CM  =  (Btu/°F)  thermal  storage  capacitance  of  body  tissues 
c  c 


CM  =  (Btu/°F)  thermal  storage  capacitance  of  skin  tissues 
s  s 


R^  -  (°F  hr/Btu)  variable  thermal  resistance  between 


conv 


ev 


rad 


w 


s  ,a 


body  core  and  skin  (see  Fig.  15) 

=  (°F  hr/Btu)  convective  heat  transfer  resistance 
between  skin  and  surrounding  air 
=  (lbwv/lbda  hr/Btu)  convective  mass  transfer 
resistance  between  skin  and  surrounding  air 
=  (°F  hr/Btu)  radiant  heat  transfer  resistance 
between  skin  and  surrounding  air 
=  (Lbwv/lbda)  specific  humidity  of  skin  and  air, 
respectively. 
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Referring  to  Fig.  16,  the  metabolic  output  (MR)  appears  as 
an  energy  source  within  the  core  of  the  body.  It  can  either  flow 
out  through  the  skin  or  be  stored  in  the  body  resulting  in  a  rise 
in  body  temperatures.  The  resistance,  R^,  is  shown  as  variable 
as  it  includes  the  effect  of  vasodilation  of  the  peripheral 
blood  vessels  with  skin  temperature  increase. 

The  thermal  capacitors  represent  the  body  energy  storage 
components.  The  core  capacitance  represents  97%  of  the  storage 
capability  and  the  skin  tissues  are  assigned  3%.  These  percent¬ 
ages  are  admittedly  assigned  somewhat  arbitrarily  for  want  of 
complete  understanding.  The  three  resistors  between  the  skin 
and  the  environment  represent  the  heat  flow  impedances  due  to 
evaporation,  convective  heat  transfer,  and  radiation.  The  radia¬ 
tion  resistance  is  shown  as  variable  to  picturize  the  transient 


thermal  radiation  pulse. 

Derivation  of  Thermal  Response  Equations 

By  consideration  of  energy  rates  at  the  nodel  points,  Tc 
and  T  ,  the  following  equations  result. 

O 

For  the  node  at  T 

c 

dTc  <TC  -  T  ) 

MR  -  C  M  -rzr- - — —  =  0  I 

c  c  dt  R 

For  the  node  at  T 

s 

T  -  T„  dT  T  -  T  w„  -  w„  T  -  T 

c  s  o»/  s  r  s  a  x  s  a  s  r  n  / 

- B - C  M  -r—  L— b -  +  K -  +  B - j  =  0  1 

d  conv  Rev  *rad 


These  can  be  rearranged 

dT_ 


cs  Ms  a 7 


dT  (T  -  T) 

cc“c3T=  MR - - 

Tc  ~  Ts  cTs  ~  Ta  +  Ms  -  wa  t  Ts  -  Tr, 
Rb  '  R=onv  Rev  _E^d 


(28) 
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Method  of  Solution 


Working  with  Eqs.  (28)  and  (29)  a  digital  computer  solution 
was  easily  obtained  bv  using  a  Runge-Kutta  method  for  solving  two 
differential  equations  simultaneously.  The  independent  variable 
was  time,  x.  The  dependent  variables  were  the  skin  temperature, 

T  ,  and  the  core  temperature,  T  .  The  initial  conditions  were: 
x  =  0,  T  =  9  5.  3° F,  T  =  9  8.1°F.  Parametric  values  were: 

w  C. 


MR 


400  Btu/hr;  body  surface  area  =  20  ft  *,  body  mass  =  160  lbm; 


body  thermal  capacitance  =  0.83  Btu/lbm0?-  h  and  h„  were  natural 

c  m 

motion  limit  values;  hc  =  .76  Btu/hr  ft^°F,  h  -  1.4  lbm  wv/hr  ft^ 
(lbm  wv/lbm  da);  T  =  92° F,  w  =  0.0129  lbm  wv/lbm  da  (note:  this 

cl  3. 

is  an  82° F  ET). 

It  should  be  noted  that  the  initial  skin  temperature  is  such 
that  the  occupant  is  actively  sweating  at  the  onset  of  the  radiant 
pulse. 

Interpretation  of  Results 

The  results  obtained  from  this  analytical  model  are  shown 

in  Fig.  17.  The  radiation  fluxes  were  taken  as  percentages  of 

curve  70-6,  Fig.  5,  Ref.  [25]. 

Radiation  -  1  =  25%  of  the  heat  flux  from  the  ceiling  is 
(see  Fig.  17)  intercepted  as  heat  flux  by  the  occupant. 

Radiation  -  2  -  double  Radiation  -  1  (i.e.,  50%  of  the 

flux  is  intercepted). 

Radiation  -  3  =  double  Radiation  -  2  (i.e.,  100%  of  the 

flux  is  intercepted) . 

Consider  the  graphical  results  for  Radiation  -  1.  It  is 
evident  that  the  occupant  can  maintain  thermoregulatory  control 
for  approximately  eight  hours  after  the  onset  of  the  radiant 
pulse  (see  skin  temperature  -  1  and  core  temperature  -  1  of 
Fig.  17).  After  eight  hours  there  is  a  loss  of  thermoregulatorv 
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control  but  the  skin  and  core  temperatures  do  not  exceed  the 
survival  limit  which  is  taken  as  a  body  core  temperature  of  102 °F. 

Referring  to  Radiation  -  2  results  it  is  seen  that  the 
person  loses  control  of  his  core  and  skin  temperature  after  about 
five  hours.  The  person  would  definitely  not  make  it  through  the 
fire  and  would  reach  a  point  beyond  the  core  temperature  limit 
in  eight  and  one-half  hours. 

Referring  to  Radiation  -  3  the  limit  would  be  reached  in 
less  than  six  hours. 

There  are  two  very  important  points  to  be  made  in  addition 
to  the  evidence  shown  in  Fig.  17.  First,  the  body  temperature 
response  of  the  shelter  occupant  will  be  strikingly  different  to 
the  radiation  pulse  depending  upon  whether  the  occupant  was 
initially  in  a  passive  moisture  loss  condition,  an  active  sweating 
condition  as  displayed  in  Fig.  17,  or  near  the  total  sweating 
condition  Cx  =  1).  Computer  runs  have  revealed  these  contrasts. 
They  are  not  reported  here  as  experimental  corroboration  of  the 
predictive  evidence  should  be  achieved  before  presenting  such 
details. 

The  second  point  has  also  been  seen  from  this  computer  model 
although  it  is  not  shown  on  the  graph.  That  is,  if  the  humidity 
is  lowered  while  keeping  the  same  effective  temperature,  the 
response  of  the  skin  and  core  will  be  relatively  slower  from  the 
same  radiation  pulse.  This  can  be  seen  by  referring  to  Eq.  (29). 
When  wa  is  lower  the  effect  of  changing  the  skin  wetness  will  have 
a  greater  effect  on  controlling  the  skin  temperature  by  evapora¬ 
tion.  Since  the  skin  wetness  is  directly  proportional  to  the 
skin  temperature  in  the  active  sweat  region,  this  will  slow  down 
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the  skin  temper?  *_re  rise  compared  to  the  higher  humidity  point 
at  the  same  effective  temperature.  This  result  is  inconsistent 
with  the  model  by  Machle  and  Hatch  which  predicts  that  the 
lowered  humidity  should  have  no  effect  on  the  time  it  takes  to 
reach  equilibrium. 

A  further  note  of  emphasis  on  the  response  of  the  body  in 
the  active  sweat  regior.  is  that  the  cooling  power  of  the  sweat 
glands  keeps  the  skin  and  core  temperatures  almost  constant 
over  a  wide  range  of  conditions.  If  the  body  is  already  heated 
to  a  point  where  the  subject  actively  sweats  and  then  the 
environment  is  changed  to  a  new  point  which  is  still  in  the 
active  sweat  region,  he  will  reach  equilibrium  very  quickly  as 
his  sweat  glands  provide  the  primary  response.  The  skin  and 
core  temperatures  will  change  only  slightly  in  reacting  to  the 
new  environment. 
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VI.  SUMMARY,  CONCLUSION  AND  RECOMMENDATIONS 
FOR  FOLLOW-ON  ACTIVITIES 

Project  Summary  and  Conclusions 

The  program  objectives  were  to  establish  a  physical  basis 
for  the  prediction  of  environmental  states  which  would  appear  as 
equally  comfortable  or  uncomfortable  to  the  human  being.  Partic¬ 
ular  attention  was  to  be  given  to  environments  with  very  low  air 
circulation  rates.  The  results  from  the  predictive  technique 
were  then  to  be  compared  wit^  past  and  current  experimental 
evidence  obtained  from  the  monoman  calorimeter.  The  results  were 
also  to  be  compared  with  evidence  found  in  the  literature.  In 
addition  to  this  main  program,  and  predicated  on  a  reasonably 
successful  derivation  of  a  predictive  model,  an  extension  of  the 
predictive  model  was  to  be  attempted.  This  extension  would 
consider  the  effect  on  the  human  subject  of  a  transient  (i.e., 
several  hours)  radiation  flux.  The  flux  would  simulate  the 
effects  of  a  building  fire  on  the  floor  above  a  shelter  space. 

The  effect  on  the  test  subject  would  ta  predicted  in  terms  of  skin 
and  body  core  time -temperature  relationships. 

The  following  summary  relates  to  these  program  objectives 
and  is  generated  from  the  previous  sections  of  this  report: 

(1)  The  relationship  of  the  percentage  of  the  body's  sweat 
glands  that  are  activated  (x)  to  the  body's  skin  temperature 
(Tg)  has  been  established  from  experimental  evidence  for  seden¬ 
tary  persons.  It  has  three  distinct  regions:  the  passive 
moisture  loss  region,  the  active  sweating  region,  and  the  region 
wherein  the  thermoregulatory  function  is  lost.  This  relation¬ 
ship  is  the  resultant  of  both  the  internal  bodv  thermoregulatorv 
processes  and  the  external  stress  caused  by  the  environment. 
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The  (x  -  Tg)  relation  is  the  physical  index  that  is  identifiable 
at  the  skin  surface  which  defines  thermal  comfort  for  sedentary 
subjects.  Values  on  this  (x  -  T  )  curve  are  identified  as 

b 

Physiological  Thermal  Index  values.  Thermal  comfort  zones  are 
also  identified  on  this  relationship.  Two  notes  of  caution; 
first,  the  need  must  be  recognized  for  more  extensive  testing 
of  subjects  especially  in  the  region  of  lost  thermoregulatory 
function  and  in  the  transitions  between  the  three  regions. 
Second,  the  relationship  applies  only  to  sedentary  persons. 
Whether  and  how  much  it  will  be  modified  for  other  activity 
levels  is  also  subject  to  further  experimentation. 

(2)  From  the  energy  balance  and  transport  rate  equations 
coupling  the  subject's  skin  surface  state  and  the  surrounding 
environment,  and  incorporating  the  (x  -  T  )  relationship,  air 
states  of  equal  comfort  (PTI  constant  lines)  can  be  established 
on  the  psychrometric  chart.  A  PTI  constant  line  will  appear  at 
different  locations  on  the  chart  depending  upon  air  flow  rates, 
metabolic  rate  per  unit  of  body  area,  amount  of  clothing,  and 
the  presence  of  thermal  radiation.  The  relocation  of  these 
lines  is  easily  accomplished  by  variation  of  computer  program 
parameters.  Hence,  comfort  zones  and  PTI  constant  lines  can  be 
established  on  the  psychrometric  chart  for  any  mix  of  these 
conditions. 

(3)  The  PTI  constant  lines  are  in  reasonable  agreement  with 
the  effective  temperature  lines  established  by  Yaglou  in  the 
82  to  85°F  ET  range.  At  lower  ET  values,  the  ET  lines  over¬ 
estimate  humidity  effects  in  comparison  with  the  PTI  constant 


lines.  Experimental  evidence  from  the  monoman  calorimeter  and 
the  work  of  other  investigators  support  the  view  that  the  PTI 
lines  are  a  better  representation  of  constant  comfort  states 
than  are  the  ET  lines  in  this  region.  For  ET  values  greater 
than  85°F  there  is  agreement  with  other  experimentors  that  in 
this  region  the  ET  lines  are  an  underestimation  of  the  humidity 
effect.  The  PTI  lines  again  agree  better  with  monoman  calorim¬ 
eter  evidence  and  with  the  other  experimental  evidence. 

(4)  The  modified  predictive  model  which  incorporates  the 
(x  -  T  )  relationship  and  body  thermal  capacitance  and  resist- 

o 

ance  effects  appears  useful  in  handling  transient  thermal  radia¬ 
tion  pulses  that  are  associated  with  fires  adjacent  to  shelter 
areas.  A  distinctly  different  body  time -temperature  response  is 
observed  depending  upon  which  of  the  three  regions  of  the 
(x  -  T  )  relationship  are  involved  during  the  transient  response. 
It  is  evident  from  these  results  that  countermeasures  for  dealing 
with  these  thermal  pulses  must  take  into  account  the  initial 
thermal  environmental  conditions  for  the  human’s  tolerance  for 
the  same  radiation  pulse  will  vary  seriously  depending  upon 
these  initial  conditions. 

The  results  generated  from  the  transient  predictive  model 
have  not  been  corroborated  by  experimental  test  work  in  the 
calorimeter. 

Recommendations  for  Follow-on  Activities 

There  are  several  areas  of  meaningful  extension  of  the 
existing  program  activities.  They  can  best  be  indicated  by 
asking  the  following  questions: 
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(1)  What  are  the  true  natures  of  the  transitional  segments  of 
the  (x  -  T  )  relationsh in  between  its  three  distinct 

D 

regions? 

(2)  Are  the  natural  motion  limiting  values  of  the  convective 
heat  and  mass  transfer  coefficients  fully  established? 

(3)  Does  the  (x  -  Tg)  relationship)  pertain  to  activity  levels 
other  than  sedentary  and  if  not  how  does  it  modulate? 

(4)  Does  a  human  respond  as  predicted  when  exposed  to 
transient  thermal  radiation  fluxes? 

There  are  many  other  questions  that  can  be  raised.  How¬ 
ever,  it  is  believed  by  the  authors  that  the  aforementioned  ones 
are  the  most  meaningful  that  can  be  answered  for  the  least  cost 
and  modification  of  the  existing  proiect  facilities. 
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In  the  earlier  stages  of  experimental  work  on  this  project, 
a  two  point  skin  temperature  harness  was  developed.  Elastic 
metallic  bands  were  used  to  hold  thermocouples  in  contact  with 
the  skin  tissues  yet  allow  evaporative  and  convective  heat 
transfer.  Because  the  early  tests  were  at  high  effective  tem¬ 
peratures,  the  skin  tissue  temperatures  were  nearly  uniform  and 
it  was  subsequently  shown  [26]  that  the  two  point  harness  pro¬ 
vided  a  satisfactory  average  value  of  the  skin  temperature,  Tg. 

Since  that  time,  environmental  test  conditions  have  been 
cooler-  and  less  humid  which  resulted  in  greater  skin  temperature 
variation  and  a  requirement  for  more  measuring  points.  The 
decision  was  made  to  use  the  same  system  of  attachment  and 
measurement  as  before  and  expand  the  number  of  skin  measuring 
points  to  twelve.  The  .twelve  points  were  the  forehead,  cheek, 
upper  arm,  forearm,  hand,  thigh,  calf,  foot,  and  four  points  on 
the  trunk.  The  average  skin  temperature  is  found  by  multiplying 
the  temperature  at  each  particular  point  by  the  percentage  of 
body  surface  area  that  particular  point  represents  (see  Ref. [27] 
for  surface  area  distribution).  For  instance,  the  calf  tem¬ 
perature  reading  would  be  multiplied  by  0.13  since  the  calf 
accounts  for  13%  of  the  skin's  surface  area.  After  a  careful 


evaluation  study  [26],  the  twelve  point  harness  was  used  for 
the  majority  of  the  tests  herein  reported. 

Problems  with  this  expanded  system  began  showing  up  almost 
.  mediately.  The  number  of  wires,  straps,  and  bands  clut  ered 
the  desk  at  which  the  test  subject  sits?  and  by  their  ve.: 


awkwardness  and  limited  freedom  precluded  any  attempts  to  prevent 
their  tangling.  The  thermocouple  wires  broke  frequently  as  a 
result  and  the  test  subjects  complained  of  the  limited  mobility 
[16]. 

As  a  result,  the  Mechanical  Engineering  Department  purchased 
a  telemetry  device — transmitter,  stepping  switch  and  batteries 
which  is  small  enough  (2”  x  1"  x  3")  to  be  worn  as  a  waist  pack, 
and  a  receiver.  The  sensors  are  thermistors.  The  thermistors 
were  incorporated  in  nylon  net  bands  that  used  velcro  tape  for 
size  adjustment  and  attachment  to  the  sites  on  the  test  subject's 
body.  This  was  a  considerable  improvement  over  the  metallic 
elastic  band  type  of  attachment  previously  used. 

The  receiver  outside  the  calorimeter  easily  monitors  the 
transmitted  signal  and  displays  the  skin  temperature  readings  t  > 
the  nearest  °C.  A  frequency  counter  permits  determination  of  the 
temperature  to  the  nearest  0.1°C.  The  stepping  switch  steps  the 
transmitted  signal  every  30  seconds.  There  are  nine  points 
incorporated  in  this  skin  temperature  measuring  harness. 

Thus,  T  is  taken  in  our  experiments  using  this  system  as 
s 

the  weighted  average  of  nine  observations  at  different  points  on 
the  body  surface  and  calculated  as  follows  [16]: 


— 

T  =  .19 

s 

Tthi  *  -14 

T  _  +  .035 

arm 

Tchk  * 

.035 

Thead 

+  .05 

Thand  *  -07 

TP  .  +  .13 
feet 

Tolf  + 

.175 

Tstch 

- 

+  .175 

Tbac!: 

where  the  readings  of  skin  temperatures  are  time  averaged  for  the 
different  parts  of  the  body. 
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APPENDIX  2 

LOGARITHMIC  MEAN  AVERAGE 

In  the  monoman  calorimeter,  the  air  temperature  and  humid¬ 
ity  changes  as  the  air  flows  along  the  human  body  length,  see 
Eig.  1. 


f  -  f 
s  I 


1  - - - - ;>.  Length 

FIG.  1:  PROPERTY  PROFILE  OF  AIR  AS  IT  PASSES  ALONG  THE 
BODY  LENGTH  AT  CONSTANT  PROPERTY  f 

5 

For  simplicity,  it  is  considered  that  under  equilibrium  condi¬ 
tions,  skin  temperature  of  the  human  body  remains  uniform  and 

constant  at  T  .  The  skin  surface  heat  and  mass  transfer  rate 
s 

equations  are  of  the  form 


where 


P  =  h  A  (f  -  f) 
p  tr  s 

| 

P  =  rate  of  transfer  of  quantity  P  per  unit  time 

A,  =  surface  area  from  which  the  transfer  occurs 
tr 

f  =  property  equilibrium  index 
f  -f  =  property  difference  which  is  the  driving  force 

3  * 

causing  transfer  of  P,  property  at  skin  surface 
minus  average  property  of  the  surrounding 
environment 

hp  =  average  convective  transrer  coefficient  for  P 
per  unit  time,  unit  area,  and  unit  property 
difference  causing  the  transfer 
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By  comparing  Eqs.  (1)  and  (6) 


<fs  -  fa>  = 


In 


fo  -  fi 

V-  'fi 

r  7 


or 


f  =  f  - 
a  s 


In  ( 


fo  -  fi 

rrrTi 


(7) 


) 


where  f  is  the  logarithmic  mean  property  value. 

cl 

In  the  case  of  convective  heat  transfer,  the  property  f  is  the 

Cl 

temperature  and  Eq.  (7)  becomes 


Ta  =  Ts  - 


T  -  T. 
o  1 


In  ( 


T  -  T. 


(8) 


±) 


Tc  -  ^ 

s  o 


and  in  the  case  of  mass  transfer,  the  property  f  is  the  specific 
humidity;  w,  and  Eq.  (7),  become 


wa  8  ws  " 


In  ( 


wo  -  wi 
v;s  -  wi 
ws  ”  wo 


(9) 


. - . . .  ... 
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APPENDIX  3 


This  program  computes  average  air  states,  sensible  and 
evaporative  heat  outputs,  inlet  and  exit  air  states  that  cor¬ 


respond  to  constant  skin  temperature  and  constant  wetness  under 
different  air  flow  rates.  The  (x  -  Tg)  relation,  shown  in 
Fig.  3  of  the  body  of  this  thesis,  is  simplified  and  represented 
by  three  straight  line  equations  representing  the  three  distinct 
regions  of  the  (x  -  T  )  curve.  The  three  straight  line  equations 
are  incorporated  in  the  program.  The  metabolic  rate  and  body 
surface  area  may  be  varied  as  program  parameters  but  are  taken 


here  for  the  average  adult  male. 
Symbols  Used  in  the  Program* 


LH  evaporative  heat  output  rate,  Btu/hr 

WI  air  humidity  ratio  at  inlet  state,  lbs  water  vapor/lb 
dry  air 

HM  convective  mass  transfer  coefficient  (determined  in  Chapter 
II,  lbm  water  vapor/hr  ft^  (lbm  water  vapor/lbm  dry  air) 

HFG  latent  heat  of  vaporization  of  one  pound  of  water,  Btu/hr 

2 

A  body  surface  area,  ft 

PVS  water  vapor  pressure  at  skin  temperature,  lbs  per  square 
feet 

WS  humidity  ratio  that  corresponds  to  saturated  air  at  T  , 
lbs  water  vapor/ lb  dry  air  c 

CP  specific  heat  of  dry  air  at  normal  temperatures, 

Btu/( lb) ( °F) 

HC  convective  heat  transfer  coefficient,  Btu/hr/sq  ft/°F 

G  mass  flow  rate,  lbm/ hr 

WO  air  humidity  ratio  at  exit  state,  lbs  water  vapor/lb  dry  air 

SH  sensible  heat  output  rate,  Btu/hr 

TI  inlet  air  dry  bulb  temperature,  °F 

TO  exit  air  dry  bulb  temperature,  °F 

TA  logarithmic  mean  average  of  air  dry  bulb  temperature,  °F 


itself. 


Symbols  not  listed  above  are  defined  within  the  program 
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WA  logarithmic  mean  average  of  air  humidity  ratio,  lbm  water 
vapor/ lbm  dry  air 

RW  water  vapor  gas  constant,  ft  lbf /( lbm) (°F) 

RA  dry  air  gas  constant,  ft  lbf /( ]bm) C°F) 

o 

P  atmospheric  pressure,  lbf /ft 
QI  air  flow  rate  at  inlet  state,  ft^/hr 

X  fraction  of  skin  surface  that  is  wet 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

1 

I 
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Program  Statements 

SUBROUTINE  GWILP  (X,TS) 

REAL  LH 

DO  13  1=1,301,25 
WI=I/10000. 

HMsl.4 

C  LATENT  HEAT  OF  VAPORIZATION  IS  TAKEN  AT  SATURATED  WATER 
TEMPERATURE  OF  95°F 
HFG=1G39 
A=  2  0 . 0 

PVS=3. 61*TS-225. 

WS=(. 622*PVS)/(2116.-PVS) 

CP=0. 24 
HC= . 765 
G=  51. 4 
Z=HM*A*X/G 

WO=WS-(WS-WI ) *EXP(-Z) 

B=AL0G( (WS-WI) /(WS-WO) ) 

LH=HM*A*X*HFG*(WO-WI) /B 
SH=  400 .  -LH 
Y=(HC*A)/(G*CP) 

TI=TS-Sf  '(G*CP*(1.-EXP(-Y))) 

TO=TI+SH/(G*CP) 

E=AL0G( (TS-TI) /(TS-TO) ) 

TA=TS-(TO-TI)/E 
WA=WS-(WO-WI)/B 
D015IG=25 ,85,15 
G—  XG 

Y=(HC*A)/(G*CP) 

Z=HM*A*X/G 

TI=TS-SH/(G*CP*( 1. -EXP(-Y) ) ) 

T0=TI+SH/ (G*CP) 

WI=WS-LH/(G*HFG*( l.-EXP(-Z) ' ) 

WO=WI+LH/(G*HFG) 

RW=1544./18.02 
RA=1544./29.95 
P=14. 7*144. 

QI=( (TI+460 . ) *G/P)  ( (WI *RW+RA) /(l.+WI) ) 

WRITE ( 5,9)X,TS,TI ,T0,WI ,WO,TA,WA,SH ,LH ,QI  ,G 
9  FORMAT  (12  F10.4) 

15  CONTINUE 

WRITE (5,11) 


mmmm 


.■*^p>fc' *2f0S?$§l 
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11  FORMAT  (1) 

13  CONTINUE 
WRITE ( 5,14) 

14  FORMAT  (11) 

RETURN 

END 

C  THE  FOLLOWING  IS  THE  MAIN  PROGRAM 

WRITE (5, 12) 

12  FORMAT( 9X , ' X ' , 8X, ' TI ’ ,  8X , 'TO* , 8X , *WI ' ,8X, ’WO* , 8X , 'TA' , 8X , 

'WA* ,8X,'SH’  8X,'LH' ,8X,'QI’) 

DO  10  ITS=  85 ,y 5 

TS=ITS 

X=( .00234+.00271*TS)/1.4 
CALL  GWILP(X,TS) 

10  CONTINUE 

DO  16  IHMX=3 ,14 
X=( IHMX*. 1) /l. 4 
TS=9  4.  87355+.  644 31* (IHMX*. 1) 

CALL  GWILP(X,TS) 

16  CONTINUE 

DO  17  JTS=96, 100,1 

TS= JTS 

X=1 

CALL  GWILP(X,TS) 

17  CONTINUE 
CALL  EXIT 
END 
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